CM multiband w7jmp 2021 march
CM Moxon rectangle by equation as descibed by
CM L.B. Cebik in www.cebik.com/moxon/moxgen.html
CM Comparison with model created using MoxGen
CM (www.qsl.net/ac6la/) yields equal results.

CE
SY Fr=14.1094			' Specify frequency
SY Fr1=21.18224			' Specify frequency 2
SY Fr2=28.21976			' Specify frequency 3
SY Fr3=43.45383			' Specify frequency 3
SY WD=1.63			' Specify wire diameter in mm/inch/WL
SY hgh=10.668			' Specify height above ground in meters (feet)

SY Q=9.999e-3	' Optimization x direction for second band
SY S=0.100346' Optimize width
SY V=0.1' optimize reflector 

SY Q1=0.012512	' Optimization x direction for second band
SY S1=0.126522' Optimize reflector

SY Inp=mm			' Enable this line for diameter in mm
'SY Inp=in			' Enable this line for diameter in inch
'SY Inp=299.7925/Fr	' Enable this line for diameter in wavelengths


SY Scal=1				' Use this line if output dimensions in meters.
'SY Scal=ft				' Use this line if output dimensions in feet.
SY WL=299.7925/(Scal*Fr)	' Wavelength in Meters/Feet
SY WL1=299.7925/(Scal*Fr1)	' Wavelength2 in Meters/Feet
SY WL2=299.7925/(Scal*Fr2)	' Wavelength3 in Meters/Feet
SY WL3=299.7925/(Scal*Fr3)	' Wavelength3 in Meters/Feet



SY DW=(WD*Inp*FR)/299.7925	' Get wire diameter in wavelengths
SY Dia=.4342945*LOG(DW)

SY DW1=(WD*Inp*FR1)/299.7925	' Get wire diameter2 in wavelengths
SY Dia1=.4342945*LOG(DW1)

SY DW2=(WD*Inp*FR2)/299.7925	' Get wire diameter2 in wavelengths
SY Dia2=.4342945*LOG(DW2)

SY DW3=(WD*Inp*FR3)/299.7925	' Get wire diameter2 in wavelengths
SY Dia3=.4342945*LOG(DW3)


'CONSTANTS
SY AA=-.0008571428571, AB=-.009571428571, AC=.3398571429

'DEFINE FIRST BAND
SY A=(AA*(Dia^2))+(AB*Dia)+AC
SY BA=-.002142857143, BB=-.02035714286, BC=.008285714286
SY B=(BA*(Dia^2))+(BB*Dia)+BC
SY CA=.001809523381, CB=.01780952381, CC=.05164285714
SY C=(CA*(Dia^2))+(CB*Dia)+CC
SY DA=.001, DB=.07178571429
SY D=(DA*Dia)+DB

'DEFINE SECOND BAND
SY A1=(AA*(Dia1^2))+(AB*Dia1)+AC
SY BA=-.002142857143, BB=-.02035714286, BC=.008285714286
SY B1=(BA*(Dia1^2))+(BB*Dia1)+BC
SY CA=.001809523381, CB=.01780952381, CC=.05164285714
SY C1=(CA*(Dia1^2))+(CB*Dia1)+CC
SY DA=.001, DB=.07178571429
SY D1=(DA*Dia1)+DB


'DEFINE THIRD BAND
SY A2=(AA*(Dia2^2))+(AB*Dia2)+AC
SY BA=-.002142857143, BB=-.02035714286, BC=.008285714286
SY B2=(BA*(Dia2^2))+(BB*Dia2)+BC
SY CA=.001809523381, CB=.01780952381, CC=.05164285714
SY C2=(CA*(Dia2^2))+(CB*Dia2)+CC
SY DA=.001, DB=.07178571429
SY D2=(DA*Dia2)+DB

'DEFINE Fourth BAND
SY A3=(AA*(Dia3^2))+(AB*Dia3)+AC
SY BA=-.002142857143, BB=-.02035714286, BC=.008285714286
SY B3=(BA*(Dia3^2))+(BB*Dia3)+BC
SY CA=.001809523381, CB=.01780952381, CC=.05164285714
SY C3=(CA*(Dia3^2))+(CB*Dia3)+CC
SY DA=.001, DB=.07178571429
SY D3=(DA*Dia3)+DB



SY rad=DW*WL/2				' Output: Wire radius in mtr/ft
SY rad1=DW1*WL1/2				' Output: Wire radius in mtr/ft
SY rad2=DW2*WL2/2				' Output: Wire radius in mtr/ft
SY rad3=DW3*WL3/2				' Output: Wire radius in mtr/ft

SY A=A*WL, B=B*WL, C=C*WL, D=D*WL	' Output: Dimensions in  meters/feet
SY A1=A1*WL1, B1=B1*WL1, C1=C1*WL1, D1=D1*WL1	' Output: Dimensions in  meters/feet
SY A2=A2*WL2, B2=B2*WL2, C2=C2*WL2, D2=D2*WL2	' Output: Dimensions in  meters/feet
SY A3=A3*WL3, B3=B3*WL3, C3=C3*WL3, D3=D3*WL3	' Output: Dimensions in  meters/feet

GW 1 35 C/2+B  -A/2 hgh C/2+B  A/2  hgh rad
GW 2 5  C/2    -A/2 hgh C/2+B  -A/2 hgh rad	
GW 3 5  C/2+B  A/2  hgh C/2    A/2  hgh rad	
GW 4 35 -C/2-D -A/2 hgh -C/2-D A/2  hgh rad
GW 5 7  -C/2   -A/2 hgh -C/2-D -A/2 hgh rad	
GW 6 7  -C/2-D A/2  hgh -C/2   A/2  hgh rad 

'Generate second band
GW 7 35 C1/2+B1  -A1/2-S hgh C1/2+B1  A1/2+S  hgh rad1
GW 8 5  C1/2    -A1/2-S hgh C1/2+B1  -A1/2-S hgh rad1	
GW 9 5  C1/2+B1  A1/2+S  hgh C1/2    A1/2+S  hgh rad1	
GW 10 35 -C1/2-D1+Q -A1/2-S hgh -C1/2-D1+Q A1/2+S  hgh rad1
GW 11 7  -C1/2+Q   -A1/2-S hgh -C1/2-D1+Q -A1/2-S hgh rad1	
GW 12 7  -C1/2-D1+Q A1/2+S  hgh -C1/2+Q   A1/2+S  hgh rad1	

'Generate third band
GW 13 35 C2/2+B2  -A2/2-S1 hgh C2/2+B2  A2/2+S1  hgh rad1
GW 14 5  C2/2    -A2/2-S1 hgh C2/2+B2  -A2/2-S1 hgh rad1	
GW 15 5  C2/2+B2  A2/2+S1  hgh C2/2    A2/2+S1  hgh rad1	
GW 16 35 -C2/2-D2+Q1 -A2/2-S1 hgh -C2/2-D2+Q1 A2/2+S1  hgh rad1
GW 17 7  -C2/2+Q1   -A2/2-S1 hgh -C2/2-D2+Q1 -A2/2-S1 hgh rad1	
GW 18 7  -C2/2-D2+Q1 A2/2+S1  hgh -C2/2+Q1   A2/2+S1  hgh rad1	

'Generate fourth band
GW 19 35 C3/2+B3 -A3/2 hgh C2/2+B3  A3/2  hgh rad1
GW 20 5  C3/2    -A3/2 hgh C3/2+B3 -A3/2  hgh rad1	
GW 21 5  C3/2+B3  A3/2 hgh C3/2     A3/2  hgh rad1	
GW 22 35 -C3/2-D3 -A3/2 hgh -C3/2-D3 A3/2  hgh rad1
GW 23 7  -C3/2   -A3/2 hgh -C3/2-D3 -A3/2 hgh rad1	
GW 24 7  -C3/2-D3 A3/2 hgh -C3/2    A3/2  hgh rad1	


GE

EK 0
TL 1 18 7 18 50 0
TL 7 18 13 18 50 0
TL 13 18 19 18 50 0
FR 0 1 0 0 Fr 0
GN 2 0 0 0 13 .005	' Average ground
EX 0 1 18 1 0
EN